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AN  APPLICATION  OF  SCMMERFELD’S  COMPLEX  ORDER  WAVE  FUNCTIONS 

TO  ANTENNA  THEORY 

Abstract 

In  the  past  wave  functions  of  integral  order  have  been  used 
quite  advantageously  ir.  the  solution  of  certain  antenna  and  boundary-value 
problems.  However,  in  some  instances  these  wave  functions  are  completely 
alien  to  the  problem  and  introduce  difficulties  which,  indeed,  can  be  re- 
solved but  only  at  the  expense  of  logical  simplicity.  To  place  in  evidence 
the  usefulness  and  "naturalness"  of  complex  order  wave  functions  for  the 
solution  of  certain  problems,  we  examine  theoretically  the  input  admittance 
of  a boss  antenna  with  the  aid  of  these  functions. 

Introduction 

Suppose  we  desire  to  construct  a solution  to  the  wave  equation 

(V2  ♦ k2)  u « 0 (1) 

where  u must  satisfy  the  Sommerfeld  radiation  condition  and  must  assume 
prescribed  values  on  the  surface  of  a sphere  of  radius  a . Using  spherical 
coordinates  r,  9,  0 with  origin  at  center  of  sphere,  assuming  that  u is 
independent  of  the  azimuthal  angle  0 , and  requiring  that  on  the  surface  of 
sphere 

a u ♦ p ^ u - f(0)  (2) 

where  a and  p are  constants,  we  consider  a solution  to  (1)  of  the  form 

oo 

v-* 

u (r.  9)  * . . A P fees  9'  h f-,1 

STo  ” » ” 

where  Pn(cos  9)  is  the  Legendre  polynomial  of  Integral  order  n and 
hn^(kr)  is  the  spherical  Hankel  function  cf  the  same  order.  The  An’s  are 
constants  which  have  to  be  determined  by  requiring  that  (3)  satisfy  (2)  when 
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r * a.  As  it  stands,  (3)  satisfies  the  radiation  condition, 

lim  r(~~  u - iku)  = 0 (U) 

r >.oc> 

because  as  r — the  spherical  Hankei  functiona  have  the  behaviour  of  an 


outwardly  propagating  wave, 


(1)  n+1  eikr 

V <*r)~<-i)  V 


To  determine  the  unknown  constants,  the  A^'SjWe  substitute  (3)  into  (2)  and 

use  the  orthogonality  property  of  the  Legendre  polynomials, 
t r 

I P„(cos  9)  P (cos  9)  sin  9 d9  • 6 — - — 

J nv  mv  nm  2m+l  (°J 

o 

where  6 ■ 0,  1 when  n / m,  n ■ m respectively.  Thus 

run 

A|n  a».i  / r(»)  P.(C03  9)  sin  » d9 

2 [a  ♦ P sir  h,n(1)(ka)]  (7l 


Substituting  (?)  into  (3)  we  have  the  desired  solution.  It  is  clear  that 
the  crux  of  this  problem  is  the  orthogonality  of  the  Legendre  polynomials. 
If  this  orthogonality  did  not  exist,  it  would  be  impossible  to  determine 
the  A^s.  However,  for  certain  problems  it  would  be  much  more  convenient 
and  "natural^  if  the  radial  functions,  instead  of  the  angular  functions, 
were  orthogonal.  This  is  trve,  for  example,  if  the  wave  function  has  to 
satisfy  on  a conical  surface  boundary  conditions  which  depend  on  r . In 


the  boss  antenna  problem  the  wave  function  must  satisfy  just  such  boundary 
conditions. 

In  this  note  we  first  develop  some  properties  of  the  complex 
order  wave  functions.  We  find  that  u can  be  represented  in  the  form  of  an 
infinite  series,  each  term  of  which  is  the  product  of  an  angular  function 
and  a radial  function.  The  angular  and  radial  functions  prove  to  be  the 
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Legendre  functions  and  the  spherical  Hankel  functions  of  complex  order 
respectively.  And  now  it  is  the  spherical  Hankel  functions  of  complex 
order  that  are  orthogonal.  Then  we  apply  these  wave  functions  to  the 
calculation  of  the  input  impedance  of  a boss  antenna. 


field, 

and 

r 

field. 

vectors 


Complex  Order  Wave  Functions 
Let  us  consider  in  spherical  coordinates  r,  9,  0,  a TM 
E ■ r Er  ♦ 9 E^  H « 0 where  r,  9,  are  the  unit  vectors 
, E0,  E.v  are  the  non-vanishing  components  of  the  electromagnetic 
For  steady-state,  i.e.,  time  dependence  exp(-icot),  the  complex 
E and  H satisfy  the  Haxwell  equations. 


V x E • 

i cop  H 

(8) 

V x H - 

i coE 

(9) 

from  which  it  follows  that 

V x V x H 

* K"  H 

(10) 

where  k2  ■ co2£p  and  6,  p denote  the  dielectric  constant  and  permea- 
bility of  free  space.  If  we  assume  that  E and  H are  independent  of  the 
azimuthal  angle  0,  then  in  spherical  coordinates  (10)  takes  the  form 


then  (11)  can  be  written  as 

ico<5  |g  [(^  ♦ k2)  u(r,  9)J  - 0 ► (13) 


We  choose  u(r,9)  such  that 

O?2  + k2)  u(r,9)  - 0 


(1U) 


for  all  permissible  values  of  r end  9 . Then  according  to  (9)  and  (12) 
the  electromagnetic  field  is  given  by  the  following  expressions: 


1 3 / < 

"sltnj  3 9 (°in 


8|?u> 


(15) 


Ui  (ru) 

r 39  a r 


v*e» 


(16) 

(17) 


Now  we  examine  the  solution  of  the  wave  equation  (lU)  that  satisfies  the 
radiation  condition  (U)  and  the  boundary  condition, 

— (ru)  ■ 0 wr  >n  r ■ a (18) 

or 

In  view  of  (16)  it  is  seen  that  the  condition  (18)  is  equivalent  to  requir- 
ing that  Eg  vanish  over  surface  of  sphere  of  radius  a * 

Writing  i’(r-0)  explicitly  in  separated  form,  n ■ R(r)  P(cos  6) 
(lU)  yields  two  ordinary  differential  equations 


1 

sin  9 


d_ 

d9 


(sin  9 ^|)  «■  CP 


0 


(19) 


1 

Tt 


d /_2  dR 
dF  (r  cBF^ 


♦ k2  r2 


c 


(20) 


where  C is  the  separation  constant.  If  we  choose  C 1)  where  ^ 

is  unrestricted  (a  complex  number),  the  solution  of  (19)  is  the  Legendre 
function  Py  (cos  9).  This  function  is  identical  to  the  hypergeometric 
function,  i.e., 

P, , (cos  9)  - F(-  1,  —'l05  9 ) (21) 

IX  - 


and  is  finite  for  all  9 except  9 ■ n . The  corresponding  solution  of 
(20)  is  the  spherical  Hankel  function  of  the  first  kind  hy^(kr)  and  it 
is  related  to  the  cylindrical  Hankel  function  of  the  first  kind  H ^^(kr) 


according  to  the  relation, 


hv(1V>  -Vjc  Hvfi/a  <kr>  <22> 

Hence,  we  can  write  the  solution  of  (1U)  as 
u(r,3)  Py(cos  9)  (23) 


where  the  summation  is  over  all  values  of  V determined  by  the  boundary  con- 
dition (18).  Substituting  (23)  into  (18)  we  obtain 

E kV  cTr  (a  ty^^ka^P^cos  9)  - 0 (2U) 

To  satisfy  (2U)  the  P'b  are  chosen  so  that 


Vie  denote  these  i^'s  as  etc. 

The  radial  functions  are  orthogonal  over  the  range 

x ■ ka  to  x •«>.  To  show  this  we  recall  from  (20)  that  for  any  P , say, 

v„ 

x (x  h^}(xj)  * (x2  - VB(  Vn  • 1))  ly^x)  - 0 (26) 


and  for  any  other  value  of  V . say,  l/ 

r IB 

X "fe  (x  hJ1)(x))  * (x2  + (x)  * 0 (27) 

m m 

Multiplying  (26)  by  h and  (27)  by  hy*^  and  then  integrating  the 
difference  from  x * ka  to  x - oo  # we  obtain 
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f* , i ( 1/  <•  i (1  V h, ,(!)#''«•>  v.  .■(!)'-*> 

"m  *'ir.  " "n  v "n  ‘ j J 'n 

J ka 


5 [‘%!1'(x)  1?  (x  h^1)(x))  - x h^m(1)  <*>  (x  hJ\,(1)<*>)]  <*  - 


tea 


d -j  oo 

X ® " X h^m  T^X^  354  (X  b^n"1^*^)  ' 

"*  J ka 


(28) 


where  the  last  equality  results  from  an  integration  by  parts.  The  integrated 
term  disappears  at  x = ka  and  x » oo  by  virtue  of  (25)  and  the  asymptotic 
behaviour  of  hi»^(i).  Thus 


co 

j hj/n(1)(x)  h)/m(])(x)  dx  - 6nm  Vn(ka) 


(29) 


ka 


Nj/n  (ka)  is  e normalization  factor  which  car.  bo  obtained  from  (28)  by  an  ap- 
plication of  1 'Hospital '3  rule  for  the  limit  jyn_ ► j/m  > and,  of  course,  6nm 
is  the  well-known  Kronecker  delta. 

Hence,  the  solution  of  (lh)  which  satisfies  the  boundary  condition 
(18)  and  tile  Sorxnerfeld  radiation  condition  is  given  by  a sum  over  all  P'bi 

u(r,9)  A j,  h^^(kr)  P^,  (cos  «)  (30) 


where  the  A^'s  are  as  vet  undetermined  constants.  Ones  we  are  given 
u(r,9^)  on  any  conical  surface  9 ■ 9^,  the  A ^ 's  are  determinable  by 
virtue  of  the  orthogonality  condition  (29).  Contrasting  (3)  and  (30)  we 
see  that  (3)  is  appropriate  for  boundary  conditions  on  a sphere,  whereas 
(30)  is  appropriate  for  cones. 

Substituting  (30)  into  (15)*  (16),  and  (17)  we  find  the  electro- 
magnetic field  in  terms  of  the  complex  order  wave  functions,  h^/^(kr)  I^*(cos  w) 


-7- 


H^Cr,©)  = lco£  / .A^hJ^kr)  ?v  (cos  9)  , (31) 

Eg(r»^)  0 (r  h^1>(kr))  ^ vcos  8)»  (3?) 

and  since  (ccs  9)  satisfies  (I9)  when  C + 1), 

Er(r,e)  - r^AJ/  15  V(1)(kr)  pj/(cos  e>  (33) 

An  Application  to  Antenna  Theory 

To  place  in  evidence  the  usefulness  of  these  complex  order  wave 

functions,  we  shall  use  them  to  compute  the  admittance  of  a boss  antenna.  A 

boss  antenna  consists  of  a coaxial  line  fitted  with  an  infinite  flange  and  a 

hemispherical  boss  at  the  end  of  the  inner  conductor  (fig.  1).  Except  for 

the  addition  of  a hemispherical  boss,  it  is  identical  to  the  circular  dif- 

(2) 

fraction  antenna. 

We  assume  the  antenna  is  excited  by  a 3ingie  propagating  mode 
(the  principal  mods)  in  the  coaxial  region.  This  mode  has  no  azimuthal 
variation.  The  non-vanishing  components  of  the  field  in  the  coaxial  region, 
z ■£»(),  a-S-p^rb,  are  H0,.  E^,  and  in  the  antenna  region,  r^a,  e< the 
non-vanishing  components  are  H^,  Er,  Eg.  Cylindrical  coordinates  are  used 
in  the  coaxial  region  and  SDherical  coordinates  in  the  antenna  region. 

In  the  antenna  region,  according  to  (33)»  the  r-component  of 
the  electric  field  is  given  by 

Er(r,9)  - V{V*  l)h^1)(kr)  ?y  (cos  9),  r>a,  9 • (33a) 

If  we  denote  Er(r,^)  by  ^(r),  (33a)  becomes 


0 - 

r2v  1)  h/j1)(kr)  P^(0), 

r>bf  9 - | 

(3Ua ) 

£(r) 

- 2v  Vj  h^(1)(kr)  py(o^ 

air^b,  9-7 
* 2 

(3Ub) 

since  the  component  of  the  electric,  field  tangent  to  the  perfectly  conduct- 
ing. baffle  mU 5 t disappear « Multiplying  both  sxdas  Ox  *-*y  KF^i  t 

and  integrating  from  kr  - ka  to  kr  = oo  , we  obtain  by  virtue  of  the  ortho- 
gonality relation  (29), 

co 

J*rlf{r)  hy  (kr)  d(kr) 


ka 


)/(Y+  1)  Py  (0)  Nw(ka) 


(35) 


Substituting  (35)  into  (31),  we  see  that  the  magnetic  field  at  any  point  in 


the  antenna  region  is  given  by 

oo 


C £ ST  ii^(kr)  h^  (kr*)  33  P„ 

V (r,9)  - i»£  J r'£(r>>  dtkr')^  „{i/t  ;)  tj/W  N,(k») 
ka 


(cos  9) 


(36) 


According  to  eq.  (2.3U)  of  reference  2,  the  magnetic  field  in  the 
coaxial  region  a£pj£b,  zfaO  is  given  by 

/\i2  - k2  ' 2 


■ ,.t  j<{»;  -■•■Sfgsp 


(37) 

where  I(z)  is  the  coaxial  line  current  and  ^(p)  is  the  p-component  of  the 
electric  field  across  the  aperture,  ^(p)  and  ^£(r)  of  (3Ub)  are  identical 
since  in  the  plane  of  the  baffle  the  r-coordinate  of  the  spherical  coordi- 
nate system  and  the  p-coordinate  of  the  cylindrical  coordinate  system  are 
identical.  The  eigenvalues  and  the  eigenfunctions  R (p)  31,(5  defined 
by 


(^~  h - + ^n2)  Rn(p)  ’ 0 


Bp 


(38) 


^ Rn(p)  " ° at  P " a’b 


and  the  A *s  are  roots  of  the  transcendental  equation 

n 


(39) 
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,T  ( X a)  K fX,b) 
O ’ ' 'n  ' o’”1* 


N (A_  a)  J (A  b) 
o ■ »*  o n 


(UO) 


The  orthogonality  relations. 


R(P)  Pdp-6, 


n"  m 


nm 


da) 


D 

J"wp) 


dp  - 0 


(U2) 


immediately  follow  from  (38)  and  (3S ) • 


Since  H0  must  be  continuous  across  the  aperture,  we  have 
* H0^"^(p,O)  for  a£p«£b.  Thus  from  (36)  and  (37) 


a 

00 


i <0 


c f r x V hy(kr)  hy  (kyt)  56  V0)  ,, 

£k  J r'f(r')  ♦ 1)  f • (0)  N (ka) 

ka  y V 


This  is  an  integral  equation.  V<e  shall  not  attempt  to  solve  it.  Rather, 
we  shall  use  it  to  formulate  a variational  principle  for  the  antenna  ad- 
mittance . 

Fhom  (2.26)  and (2.27)  of  reference  2 tho  coaxial  line  current 
and  voltage  are  given  by 


1(a)  - 2n  (a  ^kz  ♦ £eikz) 


V(z)  J Ep(~\p,z)  dp  (aeikz  - peikz)  log  | 


(Ulx) 


(U5) 


where  a and  {3  are  constants. 

The  characteristic  admittance  of  the  line  is  given  by 
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2n 


,0\ 


ft^V 


(1*6) 


Now  let  us  derive  the  variational  principle  for  the  admittance 


Y(0)  where 

Y(0)  - 

We  do  this  by  multiplying  (1*3)  by  p£(p)  and  integrating  from  p 
p - b,  and  finally  dividing  the  resultant  equation  by 


(U7) 

a to 


b 


Thus 


77  t(0)  * 


f(p)  Rn(p)  dpj2 


y^r7 


V 

[j°  £(p)dp] 


[j  £(p)<*p] 

p'dp'  ^(pOh^'1^  (kp1)^2  a 


P (0) 


1)  N^(ka) 


V0) 


(U8) 


It  can  be  easily  shown  that  Y(0)  is  stationary  with  respect  to  small  varia- 
tions of  ^(p)  about  the  true  <^(p)  determined  by  the  integral  equation  (1*3). 

We  choose  i as  a trial  function.  That  is,  we  let  &(p)  “ ~ 

p r o 

in  (1*8),  Ibe  to  the  orthogonality  relation  (1*2),  (1*8)  becomes 

b 


1(0) 


lo>£k^4j  h j' 1 ) 0<P ) dP/ 2 5Q  Pi /<°> 


2niMCktr\»V  VKp; 

" [log  |]22_»  1)  n 


,(ka) 


T-^r 


(1*9) 


'll" 


This  is  an  accurate  expression  for  the  admittance  of  the  antenna*  It  is 

sometimes  more  useful  to  consider  the  ratio  Y(0)/Y  , which  is  a dimension* 

o 

less  quantity* 


Y(0) 


> rviy» 

i (J»-y(1W>  «(■*))*  x 

^[hy^^kp)]2  d (kp) 


(|)  n (1) 

JTa" 


a p± (0> 


l)  p (0) 


/ \ 


This  is  only  a formal  solution,  however.  Since  the  values  of  the  functions 

.(3) 


involved  have  not  been  tabulated,  only  certain  limiting  values 
present  obtained  for  (50). 


can  be  at 
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Figure  1,  BOSS  ANTENNA  CONSISTING  OF  COAXIAL  LINE 
FITTED  WITH  FLANGE  AND  HEMISPHERICAL  BOSS 
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